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ABSTRACT 

We study the 7-ray emission from the pulsar magnetosphere based on outer gap mod- 
els, and the TeV radiation from pulsar wind nebulae (PWNe) through inverse Compton 
scattering using a one-zone model. We showed previously that GeV radiation fr'om the 
magnetosphere of mature pulsars with ages of ~ 10^ — 10^ years can contribute to the 
high latitude unidentified EGRET sources. We carry out Monte Carlo simulations of 
7-ray pulsars in the Galaxy and the Gould Belt, assuming values for the pulsar birth 
rate, initial position, proper motion velocity, period, and magnetic field distribution 
and evolution based on observational statistics. We select from the simulation a sam- 
ple of mature pulsars in the Galactic plane (|&| < 5°) and a sample at high latitudes 
> 5°) which could be detected by EGRET. The TeV fluxes from the pulsar wind 
nebulae of our simulated sample produced through inverse Compton scattering by 
relativistic electrons on the cosmic microwave background and synchrotron seed pho- 
tons are calculated. The predicted fluxes are consistent with the present observational 
constraints. We suggest that strong EGRET sources may be potential TeV source 
candidates for present and future ground-based TeV telescopes. 

Key words: radiation mechanisms: nonthermal - stars: statistics - stars: neutron - 
pulsars: general - 7-rays 



1 INTRODUCTION 

There are about 170 unidentified 7-ray sources in the third 
EGRET catalog, and nearly one third of these sources lie 
close to the Galactic plane |6| < 5° (Hartman et al. 1999). 
Most of those unidentified sources in the Galactic plane can 
be identified as 7-ray pulsars, possibly Geminga-like pul- 
sars which are radio quiet (Cheng & Zhang 1998; Zhang, 
Zhang & Cheng 2000). For the medium and high latitude 
sources, it has been suggested that some of them are associ- 
ated with the supernova remnants in the nearby Gould Belt 
(Gehrels et al. 2000; Grenier 2000). In addition, Harding & 
Zhang (2001) used the polar cap model (Daugherty & Hard- 
ing 1996) to investigate if 7-ray pulsars viewed at a large 
angle to the neutron star magnetic pole could contribute to 
unidentified EGRET sources in the medium latitudes associ- 
ated with the Gould Belt. Their results suggest that at least 
some of radio-quiet Gould Belt sources detected by EGRET 
could be such off-beam 7-ray pulsars. 

At the same time, these 7-ray pulsars could produce 
wind nebulae through the interactions between relativistic 
wind particles with the interstellar medium (ISM). The pul- 
sar wind nebulae will contribute to the production of non- 
pulsed X-ray emission by synchrotron processes (Chevalier 



2000), and TeV photons through inverse Compton scatter- 
ing (ICS)(Aharonian, Atoyan & Kifune 1997). These ex- 
cess TeV photons have been detected in some known pulsar 
wind nebulae, such as the Crab, Vela, PSR 1706-44, possibly 
Geminga (Kifune et al. 1995; Yoshikoshi et al. 1997; Aha- 
ronian et al. 1999; Lessard et al. 2000). Therefore, if 7-ray 
pulsars contribute to the unidentified EGRET sources, pos- 
sible TeV signals could be expected to be detected in these 
EGRET sources. Several groups have searched for TeV sig- 
nals in the error boxes of unidentified EGRET sources, for 
example, with the HEGRA AIROBICC array (Aharonian et 
al. 2002a), and the Whipple 10m Gamma-Ray Telescope. No 
TeV source detection has been confirmed at Whipple, with 
only an upper limit TeV fiux for about 20 EGRET sources 
determined at ~ (3 — 6) x 10~^^ photon cm~^ s~^ (Fegan 
& Weekes 2004). Deep observations of the Cygnus region 
using HEGRA (Aharonian et al. 2002b) showed an uniden- 
tified TeV source in the vicinity of Cygnus OB2 with an 
integral fiux F{> 1 TeV) ~ 5 x 10"" photon cm"^ s"\ at 
the edge of the 95% error circle of the EGRET source 3EG 
J2033-I-4188. However, no connection between the EGRET 
source and the TeV signal has yet been confirmed. 

It has been shown that GeV photons can be produced 
in the pulsar magnetosphere in outer gap models (Cheng et 
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al. 1986; Zhang & Cheng 1997). A revised outer gap model 
(Zhang et al. 2004) takes into account the effect of the incli- 
nation angle a between the magnetic axis and the rotational 
axis, which can determine the gap size of the outer gap. 
This allows some pulsars with appropriate combinations of 
a, P and B, to maintain the outer gap for at least ~ 10® 
years. Their advanced ages allow these pulsars enough time 
to move up to high Galactic latitudes as weak 7-ray sources. 
This leads Cheng et al. (2004a) to propose that mature 7- 
ray pulsars with ages ~ 10® — 10® years can contribute to 
the unidentified EGRET sources. These mature pulsars also 
remain active in producing relativistic wind particles, and 
form compact wind nebulae. In addition, TeV photons can 
be created in the nebulae through the ICS process. In this 
paper, we study the possible connection between TeV 7-ray 
sources and the Unidentified EGRET sources. We do not 
know if the Unidentified EGRET sources are pulsars; even 
if they axe pulsars we still do not know their properties, 
i.e. period, magnetic field, inclination angle, distance etc. 
Without these parameters we cannot calculate their 7-ray 
properties. Therefore, wc apply a statistical approach using 
Monte Carlo simulation to study the Unidentified Gamma- 
ray EGRET Sources. First we will simulate the galactic pul- 
sar population and use the outer gap model to calculate the 
MeV-GeV photon power from these simulated pulsars. We 
have ignored the contribution from the polar gap for simplic- 
ity. We can determine which simulated pulsars can be de- 
tected by EGRET in 7-rays; we call them 7-ray loud pulsars. 
The next step is to calculate the 7-ray emission from the 
pulsar wind based on the simulated pulsar parameters. We 
should point out that the distribution of 7-ray loud pulsars 
is model dependent. Subsequently we study the TeV 7-rays 
emitted from the pulsar wind when they interact with their 
ambient interstellar medium. We argue that strong EGRET 
sources may be potential TeV source candidates for current 
and future TeV telescopes. 

In §2, theories of radiation of mature pulsars from both 
inside and outside the light cylinder will be briefly reviewed. 
A detailed discussions of the ICS processes of relativistic 
electrons on the background and synchrotron photons us- 
ing a one-zone wind nebula model (Chevalier 2000; Cheng, 
Taam & Wang 2004b) will be presented in § 3. In § 4, the 
properties of the identified TeV pulsar wind nebulae sources 
arc reviewed and compared with our models. In § 5, we de- 
scribe the Monte Carlo simulations, similar to those carried 
in Cheng et al. (2004a), for deriving the distribution of 7- 
ray pulsars in the Galaxy and Gould Belt which can be 
detected by EGRET. We then calculate the expected GeV 
and TeV fiuxes from these 7-ray pulsars generated by the 
Monte Carlo simulations. The TeV flux distribution of the 
unidentified EGRET source candidates, which could possi- 
bly be detected by the present and future ground-based TeV 
telescopes, will be presented in § 6. Finally, we present our 
summary and discussions in § 7. 



2 RADIATION THEORIES ASSOCIATED 
WITH PULSARS 

In this section, we will review the emission properties of 
mature pulsars whose ages axe from ~ 10® — 10® yeaxs. The 
high energy radiation of mature pulsars can come both from 



the pulsar magnetosphere inside light cylinder, and from the 
pulsar wind nebula outside light cylinder. 

2.1 Inside light cylinder - outer gap model 

We present the 7-ray emission properties of pulsars us- 
ing the outer gap models originally proposed by Cheng 
et al. (1986a,b). Based on the model, Zhang & Cheng 
(1997) have developed a self-consistent mechanism to de- 
scribe the high energy radiation from spin-powered pul- 
sars. In their model, relativistic charged particles from a 
thick outer magnetospheric accelerator (outer gap) radiate 
through the synchro-curvature radiation mechanism (Cheng 
& Zhang 1996) rather than the synchrotron and curvature 
mechanisms in general, producing non-thermal photons from 
the primary pairs along the curved magnetic field lines 
in the outer gap. 

The characteristic emission energy of high energy pho- 
tons emitted from the outer gap is given by (Zhang & Cheng 
1997) 

E-,,, ~ 5 X l0V^/2Bf/^P~^/*(^)-«/« eV, (1) 

where P is the rotation period, S12 is the dipolar magnetic 
field in units of 10^^ G, Rl = cP/2ir is the light cylinder 
radius, r is the distance to the neutron star, and / is the 
fractional size of the outer gap. The 7-ray spectrum drops 
exponentially beyond the energy Ej^c- 

The factor /, defined as the the ratio between the mean 
vertical separation of the outer gap boundaries in the plane 
of the rotation axis and the magnetic axis to the light cylin- 
der radius, is limited by the pair production between the 
soft thermal X-rays from the neutron star surface and the 
high energy 7-ray photons emitted from the outer gap re- 
gion, and can be approximated as / ~ 5.5P^'^^^^B^2 . The 
size of / in turn determines the total 7-ray luminosity of the 
pulsar, which is given by (Zhang & Cheng 1997) 

— f^Lsd, (2) 

where Lgd = 3.8 x 10'^^Bi2i-'~* erg is the pulsar spin 
down power. However, the estimation of gap size / by Zhang 
& Cheng (1997) does not include the effect of inclination an- 
gle and their model can produce very few pulsars with age 
~ 10'' — 10" yrs. Zhang et al. (2004) have studied the prop- 
erties of the outer gap by including the effect of the inclina- 
tion angle. Also, instead of taking half of the light cylinder 
as the representation of the outer gap, they have calculated 
the average outer gap size and use it as the representation of 
the outer gap. This is very important for a statistical study. 
Cheng et al (2004a) show that in Monte Carlo simulations 
even if the inclination angle of pulsars is randomly selected 
from a uniform distribution, the simulated 7-ray pulsars de- 
tected by EGRET in the galactic plane are younger and 
tend to have a larger inclination angle. On the other hand, 
7-ray pulsars detected by EGRET at higher latitudes are 
older and have a smaller inclination angle. In other words, 
even if the seed distribution has a uniform distribution, a 
non-uniform distribution can be generated by the 7-ray se- 
lection effects. In this paper, we intend to use Monte Carlo 
methods to study GeV and TeV 7-ray properties of pulsars, 
so we adopt the model of Zhang et al. (2004) to determine 
the size of the outer gap. Assuming that the representative 
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region of the outer gap is the average distance to the gap, 
the mean fractional size of the outer gap can be approxi- 
mated as f{a,P,B) ^ Ti{a,P,B)f{P,B), whore 7j{a, P, B) 
is a monotonically increasing function of a, B, and P. The 
value of rj roughly decreases by a factor 3 from large incli- 
nation angles to smaller angles (Zhang et al. 2004). These 
7-rays from the pulsar magnetosphere will contribute to the 
pulsed GcV photons in 7-ray pulsars detected by EGRET. 

In general the differential energy spectrum of 7-rays for 
each pulsar is different. In principle we can calculate it if the 
pulsar parameters are specified. However, it is very difficult 
to do so in a Monte Carlo simulation, in which wo need 
to deal with over ten million pulsars. EGRET has detected 
six 7-ray pulsars and their energy spectra from 100 MeV 
to a few GeV are all very close to ( Hartman et al. 
1999). For thin outer gap pulsars, it has been shown that 
the energy spectrum is proportional to E^^ log{Emax / E^) oc 
E~^ (Cheng, Ho and Ruderman 1986b; Cheng and Ding 
1994). For simplicity in the Monte Carlo simulation, we will 
approximate the expected energy differential 7-ray flux of 
the pulsar as 

where d is the distance of the pulsar, and Af2 is the the 
solid angle of 7-ray beaming. The value of Afl generally 
varies with different pulsars. For simplicity, we assume a 
constant beaming solid angle Afi ~ 1 sr in all our analyses. 
In order to compare our model results with observations, we 
calculate the expected integral EGRET flux of our model 
7-ray pulsars using the formula 

r^max 

Sj{E^ > lOOMeV) = / F{E^)dE^, (4) 

J IOOMgV 

where Emax is the maximum 7-ray energy detected by 
EGRET, chosen to be 50 GeV. 

Here, we have ignored the 7-ray contribution from 
the polar cap, which could be important for Unidentified 
EGRET Sources (Gonthicr ct al. 2002). However, the phase- 
resolved EGRET data of the Crab pulsar (Cheng , Ruder- 
man and Zhang 2000), Geminga (Zhang and Cheng 2001) 
and the Vela pulsar (Romani 1996) can be explained very 
well by the outer gap model and it appears that polar cap 
emission is unimportant. Most recently, Muslimov & Hard- 
ing(2004) have suggested that the slot gap model could be 
able to explain the phase-resolved data as well. Since the 
distributions of simulated 7-ray pulsars are model depen- 
dent, in this paper we will ignore the contribution of 7-rays 
from the polar cap for simplicity. 

2.2 Outside light cylinder - pulsar wind nebulae 

Previous models of the pulsar wind nebulae (e.g.. Kennel 
& Coroniti 1984; Chevalier 2000) mainly concentrate on the 
bright nebulae produced by interactions between young pul- 
sar wind particles and the supernova remnant (SNR). On 
the other hand, while there arc no SNR surrounding the 
mature pulsars, they remain active enough to produce the 
relatively faint, compact synchrotron nebulae through inter- 
actions between the relativistic wind particles and the inter- 
stellar medium. In this work, we use the one-zone model 
to describe high energy radiation from pulsar wind nebulae 



(Chevalier 2000). In the model, the relativistic electrons in 
the shock waves emit X-rays through synchrotron radiation. 
This model can well explain the X-ray luminosity and spec- 
tral properties of pulsar wind nebulae (Cheng et al. 2004b). 

Mature pulsars move at a high proper velocity after 
their birth, and can form a bow shock structure to produce 
synchrotron wind nebulae when the pulsar proper motion 
velocity is larger than the sound speed in the ambient inter- 
stellar medium (ISM). The characteristic size of the shock 
wave produced by interactions between the pulsar wind par- 
ticles and the ISM is referred to as the termination radius, 
Ra, and can be derived from 

D fT In '2 \l/2 iril6rl/2 -1/2 -1 r-\ 

Rs = {Lsd/2npvpc) ~ 10 L^^^^^n Wp^ggo cm, (o) 

where p = nrup, n = 1 cm~^ is the number density of the 
ISM, rUp is the proton rest mass, Vp is the pulsar velocity in 
units of 350 km s~^ (cf. § 5.2), and Lsd,34 is the spin down 
power of the pulsar in units of 10^'* erg s~^. 

We have assumed the pulsar wind can carry away most 
of the pulsar spin-down power (Lgj) and deposit in the shock 
waves when the wind interacts with the interstellar medium. 
In general, the energy in the shock waves is stored in the 
magnetic field, the energetic protons (ions) and electrons. 
However, the fractional energy density of the magnetic field 
es is low (typically, es ~ 0.001 — 0.01, Kennel & Coroniti 
1984), then by assuming cquipartition of energy between 
protons and electrons, we obtain tp ^ te ^ 0.5. For a given 
es, the magnetic field at the termination radius is estimated 
as B = (6esI/sd/i?sC)^^^. At the shock front, the electron 
energy distribution is N{^') (x 7"^ for 7™ < 7 < 7max, 
where 7^ ~ [{p — 2)/(p — l)]r,'7ra, and 7u, is the Lorentz 
factor of the pulsar wind particles. 

We derive the value of 7max by two methods. First, 
an estimate for 7max can be obtained by equating the 
synchrotron cooling timescale to the electron accelera- 
tion timescale. The former timescale is given by tsyn = 
Q-KrUec/ aT^B^ , and the latter timescale is given by tacc ~ 
7mec/eB, leading to 7max,i = {Q-Ke/ (jtBY'"^ ~ lQ^°BZl''\ 
where ot is the Thompson cross section, and B-s denotes 
B/10~^G. Second, we can estimate 7max by equating the 
acceleration timescale and diffusion timescale of electrons. 
The electron diffusion time in the termination radius can be 
estimated as tdift ~ {Rs/rLfrL/c, where tl ^ 'ymeC^ /eB is 
the Larmor radius. This gives a maximum Lorentz factor of 
7max,2 ~ 10*B-5i?.,i6, where Rs,w = Rs/lO"^ cm. We will 
assume 

7max = min(7inax,l, 7max,2) (6) 

in the following discussion. 

The Lorentz factor of wind particles cannot be too 
small, the production of the synchrotron nebulae requiring, 
for example, 7„ > 10''. Theoretical constraints on 'y^ arc 
still difficult to determine at present. In this paper, we use 
the arguments of Rudermau (1981) and Arons (1983) that 
large fiuxes of protons (ions) can be extracted from neutron 
stars and can be accelerated in the parallel electric field in 
the magnetosphere. Wc also assume that the initial Poynt- 
ing flux can be converted into thermal and kinetic energy 
of the particles well within the termination radius, and that 
the Lorentz factors of electrons and protons are the same. 
For young pulsars like the Crab pulsar and the Vela pulsar, a 
large number of electrons and positrons can be produced in 
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the outer magnetosphere (Cheng, Ruderman & Ho 1986b) 
and hence electrons/positrons are the main energy carriers 
in the pulsar wind. However, for mature pulsars with age 
> 10^ years, the outer gap cannot produce too many pairs 
(Zhang & Cheng 1997). The pair creation will be mainly 
at the polar cap, however, the ratio between electron num- 
ber and proton number is < nip /me, in particular for the 
space-charge-lirnited flow situation (Arons & Scliarlemann 
1979; Cheng & Ruderman 1980). Since we will show that 
the main contributors of TeV 7-ray sources are mature pul- 
sars, so we shall assume that protons are the main energy 
carriers of the pulsar wind. We want to point out that this 
assumption does not conflict with the equipartition assump- 
tion in the shock waves because pulsar wind energy will be 
redistributed when the wind interacts with the interstellar 
medium in the shock waves. 

The protons can then carry away most of the spin down 
power in this scenario (Coroniti 1990). Thus, we can de- 
rive the spin down power as L^d ~ Njw'm.pC^ , where A'" 
is the outflow current from the surface. This outflow cur- 
rent should be of the same order as the Goldreich- Julian 
current (Goldreich & Julian 1969) which is given as iV ~ 
1.35 X 10^°Si2P"^ s"^ oc lI^"^. For mature pulsars with 
Lsd ~ 10^''"^® erg s~^, we derive atypical value of jui ~ 10^- 

The critical synchrotron frequency of relativistic elec- 
trons with a Lorentz factor 7 is 1^(7) = 7^eS/27rmeC. The 
synchrotron radiation properties from the pulsar wind neb- 
ulae depend on two characteristic frequencies: Um, radiated 
by the electrons with Lorentz factor 7™, and the cooling fre- 
quency, fc, given as fc = e/{2-!TmecB^)[{6nmeC)/{aTto)]^, 
where to is the characteristic timescale of the nebula es- 
timated from the flow timescale in a characteristic radia- 
tion region. If we define a characteristic time scale to ~ 
Rs/vp ~ l(f s, and B ~ 10~® G, then the cooling fre- 
quency Vc ~ 10^^ Hz. This corresponds to a characteris- 
tic cooling Lorentz factor 7c ~ W'' , which is significantly 
larger than 7^- Therefore, we can derive the luminosity and 
spectral properties of synchrotron radiation from the pul- 
sar wind nebulae. When Um < v < Vc, it can be described 
as a non-thermal spectrum with a power law photon index 
r ~ (p+ l)/2, and a synchrotron radiation luminosity esti- 
mated by (also see Chevalier 2000; Cheng et al. 2004b) 



J^syn^fj — — — TTTT; 



9mi''-''/''c(p-i)/2 

Vp 7u, lo-t^sdl' 



(7) 



This case is called the slow cooling regime. On the other 
hand in the fast cooling regime when v > i/c, the emission 
spectrum will steepen because of the cooling of relativistic 
electrons, the photon index F ~ (p+2) /2, and the luminosity 



L.,„(j.) ^ i(|5|)--i(-£-)(--2)/2(nmp)(-^>/\r^ 



3 TEV PHOTONS FROM PULSAR WIND 

NEBULAE THROUGH INVERSE COMPTON 
PROCESSES 

The synchrotron radiation of relativistic electrons in the pul- 
sar wind nebulae can well contribute to non-thermal radio 
to X-ray emissions. At the same time, the inverse Comp- 
ton scattering (ICS) of the same relativistic electrons on 
the ambient photon fields results in the production of the 
very high energy TcV photons observed (de Jager & Harding 
1992; Atoyan & Aharonian 1996; Aharonian, Atoyan & Ki- 
fune 1997). In this section, we will discuss two cases of the 
ICS processes of two different types of seed photons, namely, 
background and synchrotron photons. We can see that the 
Inverse Compton luminosity is sensitively dependent on the 
properties of the seed photons. 

3.1 Inverse Compton Scattering from background 
photons 

The relative magnitude between the emission from syn- 
chrotron and ICS processes depends on the magnetic field 
B and the photon field density Wph. The magnetic field 
in the interstellar medium is typically within the range 
(0.3 - 1) X 10"^ G. Synchrotron X -ray photons and 7-ray 
photons firom ICS have been shown to be produced in re- 
gions of different magnetic fields (Aharonian et al. 1997). In 
this paper, we assume for simplicity that the X-rays and 
7-rays come from the same region with a magnetic field 
B ~ 10"^ G. 

The photon energy density surrounding the pulsars can 
be attributed to several origins: first, the cosmic microwave 
background (CMB) radiation contributes a constant en- 
ergy density wcmb ~ 4 x 10~'^^erg cm~^, second, the dif- 
fuse galactic dust far-infrared (FIR) and star light near- 
infrared and optical background (si), and third, possible 
radiation fields of local origin. The density of the galactic 
background field varies from site to site, with an average 
values wfir ~ 10~^^ erg cm~^ and Wsi ~ lO"'^^ erg cm~^ 
(Mathis, Metzger & Panagia 1983). Although the energy 
density of the galactic radiation background is higher than 
that of the CMB, its importance is strongly reduced in the 
energy region above 100 GeV due to the Klein-Nishina ef- 
fect (Aharonian et al. 1997) . Therefore we consider here only 
the CMB photon field as the external seed photons for the 
production of TeV photons. 

The characteristic energy, Esyn, of synchrotron radia- 
tion from relativistic electrons with typical Lorentz factor 7 
is given as 



?syn = hvsyn = ^7 eB / nieC ~ 1_B_5(7/10 ) keV. 



(9) 



The characteristic energy of the photons from the ICS of 
these same energy electrons on the CMB photons is given 
by 

Sic = -r'hvo ~ l(7/10*)''(/ii'o/ecMB) TeV, (10) 

where ecmb ~ 10"* eV is the energy of the CMB photons. 
The ratio of the luminosity of the emission from the ICS 
process to that from synchrotron radiation can be simply 
obtained as {jhuo < nieC^) 

Lie _ Wcmb 

Lsyn WB 



(11) 
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where wb = B^/8n is the energy density of the magnetic 
field. To derive the 7-ray luminosity of TeV photons by the 
ICS process, wo can first estimate the synchrotron X-ray 
luminosity of the pulsar wind nebulae using the one-zone 
model discussed in § 2.2. 

The TeV photon luminosity and flux distributions will 
be studied in detail in § 6 using a simulated sample of mature 
pulsars. Here, we use some typical pulsar parameters for a 
simple estimate. According to equations (9) and (10), TeV 
photons through ICS produced by relativistic electrons cor- 
respond to X-rays (~ 1 keV) through synchrotron processes. 
In the calculation, f ~ 10^* Hz denotes the X-rays through 
the synchrotron process, and we take some typical param- 
eters of pulsars and nebulae to find the TeV luminosity: 
L,d ~ 10^" erg s-\ee ~ 0.5, es ~ 0.01, 7„ ~ 10^ J?^ ~ 10^® 
cm, the electron energy spectral index p — 2.2, then we 
obtain the TeV luminosity by ICS processes Lie (TeV) ~ 
10^1 erg s-^ 

The spectral profile is also related to the synchrotron 
process. As discussed in §2.2, there exist three Lorentz 
factors of the relativistic electrons 7m,7mai, and 7c. The 
Lorentz factor of the ICS electrons should satisfy yhvcMB < 
meC? , which gives a critical Lorentz factor ^cut ~ 10^. Then 
when the Lorentz factors satisfy the relation 7^ < 7c < 
7ciit < 7maa;, for ICS elcctrous with 7 between 7^ and 7c, 
the ICS spectrum will show a power-law with a photon index 
r = (p + l)/2, for 7 between 7c and 7cut, the photon index 
becomes V = {p + 2)/2. For ICS electron 7 above 7cut or 
'ymax, the spectrum will become am exponential decay form 
reflecting a high energy cut-off. For mature pulsars discussed 
above with ')max ~ 10*, a cut-off energy at ~1 TeV will be 
produced. 

3.2 Synchrotron self-Compton scattering 
contribution 

Using the synchrotron spectrum as the sources of seed pho- 
tons, we can compute the resulting ICS emission (called syn- 
chrotron self-Compton, SSC, photons). Similar to the origi- 
nal synchrotron spectrum, the up-scattered component also 
has two characteristic energies, ~ /i7^i'm,, and E° ~ 
h-y'^i'c- In addition, the condition 7/11^(7) < rrieC^ will give a 
characteristic critical factor 7crit ~ 2 x 10^ for B = 10~^ G 
and a critical energy i?crit ~ 7crit'"ieC . Due to the limitation 
that ^maxhvo < meC^ , we can find a frequency vo such that 
if 1^0 > ~ 50(7„a^/10*) TeV. In the 

case of V{) < Vrn, a lower limit of "fmhi^m ^ rrieC^ is required, 
thus implying a cut-off energy E^^^ ~ j'^nieC^ ^ 10 TeV 
for B ~ 10"' G and 7™ ~ 10^ If Efif^ < Ef^ < E^^^, 
then the SSC spectrum will show a power-law with a pho- 
ton index T = {p + l)/2 for photon energy below Sf^*^, 
and an index T — (p + 2)/2 for photon energy above E^^'~^. 
For photon energy > iJ ~ "fcritmeC^ , the energy spectrum 
of SSC becomes steeper because of the Klein-Nishina effect 
(Aharonian et al. 1997). 

In general, the intensity of the ICS emission relative 
to the synchrotron radiation can be estimated by consid- 
ering the ratio of the energy density of synchrotron ra- 
diation and the magnetic field, Wsyn/ws. In our one-zone 
model, synchrotron radiation comes from the termination 
radius Rs- We approximate the synchrotron energy density 
as WJgyn ^ 3Lsyntf/4-KRl, where tf ~ Rs/c. The SSC pho- 



ton energy is given by Esse ~ l^hv, where v ~ ■y^eB/meC. 
Because the relativistic electrons producing the TeV pho- 
tons by SSC correspond to a characteristic Lorentz factor 
7 ~ 10® < 7c (assuming B ~ 10"' G and to ~ 10^ s), 
we should apply equation (7) to compute the energy density 
and luminosity of low energy synchrotron seed photons. The 
ratio of the SSC and synchrotron luminosities can be given 
by 

Lssc ^ aTe<^-^>/"(nmp)<''+^)/^ ^p_i^(p-3)/4 
isyn 9ml''-''''^c(f-i)/2 

.r^^/^^^^oe-//^ (12) 

Assuming i^sccd ~ 10^'* Hz, to ~ 10^ s, and taking the 
other parameters the same as those in § 3.1, we determine 
the luminosity of TeV photons produced by SSC as about 
issc(TeV) ~ 10^° erg s-\ 



4 IDENTIFIED TEV SOURCES OF PULSAR 
WIND NEBULAE 

The present TeV telescopes have detected excess TeV pho- 
tons from three 7-ray pulsars, namely. Crab, Vela and PSR 
1706-44, with a possible detection from Geminga (see re- 
views by Aharonian 1999; Weekes 2004). Their 7-ray emis- 
sions in the GeV band are dominated by contributions from 
the pulsar magnetosphere, as predicted in the outer gap 
model (§ 2.1). However, no pulsed TeV photons are detected 
by the present telescopes (Kifune et al. 1995; Yoshikoshi et 
al. 1997; Aharonian et al. 1999; Lessard et al. 2000), thus 
implying that the TeV photons mainly come from their wind 
nebulae. In this section, we check our one-zone model by es- 
timating the TeV signals from these four sources. Because 
of the simplified nature of these models, we do not expect 
to reproduce the detailed properties of the nebula. More de- 
tailed discussions for individual sources have been addressed 
by other authors (e.g. de Jager & Harding 1992; Atoyan 
& Aharonian 1996; Aharonian, Atoyan & Kifune 1997). 

4.1 Crab nebula 

The Crab nebula is the first pulsar wind nebula with TeV 
signal detection (Weekes et al. 1989), and it has been ex- 
tensively studied by many TeV detectors since, for exam- 
ple, Milagro (Atkins et al. 2003), HEGRA (Aharonian et 
al. 2000), Whipple (Lessard et al. 2000), Tibet (Amenomori 
et al. 1999), and CANGAROO (Tamimori et al. 1998). De- 
tailed spectral and flux information have been obtained. The 
Crab TeV spectrum can be fitted in the energy range above 
1 TeV by a simple power law with a photon index ~ 2.5±0.2, 
and shows no cut-off up to 100 TeV. The TeV luminosity is 
estimated to be ~ lO^** erg s"^. 

Atoyan & Aharonian (1996) discussed the radiation 
mechanism of the Crab nebula based on the model of Ken- 
nel & Coroniti (1984), and fitted the observed data from 
radio to GcV by synchrotron radiation, and the TeV band 
through ICS processes. With a complicated spatial struc- 
ture, the Crab nebula probably cannot be fully described 
in our one-zone model. What we present in the following 
should therefore be treated as a rough check. 

The magnetic field in the Crab nebula is relatively high 
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at B ~ 10"^ - 10"" G (Hester et al. 1996), which signifi- 
cantly reduces the relative contributions from the ICS pro- 
cesses. With a ratio wcmb/wb < 10~*, the SSC contribu- 
tion to TeV photons should be the dominant component. 
Assuming B ~ 10~* G, we derive the cooling Lorentz factor 
of electrons to be jc ~ 10®, and the characteristic timescale 
to ~ 3 X 10^° s. The imphed cooUng energy ~ 1 TeV, 

which is only 1% of the predicted SSC cut-ofT. Therefore, 
above 1 TeV, the photon index of the ICS spectrum is given 
by r = (p -h 2)/2, where p > 2. On the other hand, the ob- 
served TeV spectrum of the Crab nebula implies that p ~ 3. 
However, from the observed synchrotron X-ray spectrum of 
the Crab nebula (Willingale et al. 2001), F = (p+2)/2 ~ 2.1, 
which gives p ~ 2.2. Since jcrit ~ 7c ~ 10® in the Crab neb- 
ula, and above ■jcritmeC^ ~ 1 TeV, the SSC spectrum can be 
steeper than (p + 2)/2, hence agreeing with the observations. 
The discrepancy may also result from the spatial variations 
of the electron spectrum and magnetic field in the Crab neb- 
ula considered by other models (de Jager & Harding 1992; 
Atoyan & Aharonian 1996; SoUerman et al. 2000), which are 
not considered in our one-zone model. 

Despite the very large size of the Crab nebula, the 
TeV photons come from a central region ~ 1 pc (Atkins 
et al. 2003). The termination radius of the nebula is about 
Rs = A X W^^ cm. In our one-zone model, we can take 
te es ~ 0-5, p — 2.2, and 7„ ~ 10^, to estimate the 
synchrotron luminosity and energy density. At ~ lO'^'* Hz, 
the luminosity Lsyn ~ 10^® erg and Wsyn/ios ~ 0.03. 
Therefore the SSC contribution to TeV luminosity is Lssc ~ 
3 X 10** erg s~^, which is comparable to the observed value. 

4.2 Vela and PSR 1706-44 

Yoshikoshi ct al. (1997) used the 3.8 m imaging Ccrcrikov 
telescope near Woomcra, South Australia to detect TeV pho- 
tons of the Vela pulsar region. The detected TeV 7-ray flux 
is 2.9 X 10~^^photon cm~^ s~^, corresponding to a lumi- 
nosity of 3x10*^ erg s~^ assuming the distance of Vela to 
be 300 pc. Similar to Vela, PSR 1706-44 has also been de- 
tected in the TeV energy range (Kifune et al. 1995; Chad- 
wick et al. 1998). The fiux above a threshold of 1 TeV is 
~ 10~^^ photon cm~^ s~^. However, recent observations of 
HESS reported non-detection of PSR 1706-44, with an up- 
per limit flux --^ 1.3 x lO^^^erg cm^^ s^^ (Aharonian et al. 
2005). The reason for the observational discrepancy is un- 
known, possibly the TeV emission could be variable on a 
time scale of years, e.g. in the different medium environ- 
ments. Assuming a distance of 1.5 kpc, the TeV luminosity 
is lower than 4 x 10^^ erg s~^. TeV spectral information has 
been obtained for neither of these two sources. 

Vela and PSR 1706-44 have nearly the same spin down 
power of a few times 10^® erg s~^. This suggests that their 
nebulae have similar properties moving in the interstellar 
medium. Assuming te ~ 0.5 and es ~ 0.001, wo estimate the 
magnetic field B in the two pulsar wind nebulae. According 
to the Chandra observation of Vela (Pavlov et al. 2001), 
together with Rs ~ 10^^ cm (Cheng et al. 2004b), and L^d ~ 
6.9 x 10*® erg s~\ we estimate that B ~ 10"^ G. For PSR 
1706-44, using Lad ~ 3.4x10*® erg s"\ and Rs ~ SxlO^'' cm 
(Finley et al. 1998; Cheng et al. 2004b), we estimate B ~ 
3 X 10~® G. The observed nebula photon indices of Vela and 
PSR 1706-44 in X-rays is about 1.7, which requires Vc > I'x, 



and p ~ 2.4. As discussed in § 3.2, the ICS processes by 
CMB photons would be the dominant contributor to the 
TeV flux of the two pulsar wind nebulae. The proper motion 
velocity of Vela is about 65 km s~^ (Pavlov et al. 2001). 
Assuming for the two pulsars to ~ 10^® s, ~ 10®, and u ~ 
3 X lO^'^ Hz, to find the synchrotron radiation, we estimate 
the TeV luminosities of the pulsar nebulae to be ~ 5 x 
10*^ erg for Vela, and ~ 2 x 10** erg s'^ for PSR 1706- 
44. So the present assumption of model parameters for two 
pulsar wind nebulae may predict higher luminosities than 
are observed. 

4.3 Geminga 

Geminga is a strong 7-ray pulsar source in the GeV energy 
range, but is very weak in other energy bands. For example, 
its X-ray luminosity (including pulsar and the wind nebula) 
is about ~ 10*® erg s"'^ (Caraveo et al. 2003). Observations 
of Geminga have not been yielded a significant excess of 
TeV photons. Aharonian et al. (1999) obtained an upper 
limit which is only ~ 13% of that of the Crab nebula, cor- 
responding to a luminosity of < 6 x 10*® erg s~^ assuming 
a distance of 160 pc. 

The wind nebula around Geminga has a typical bow 
shock structure, with a non-thermal X-ray spectrum and 
a photon index T ~ 1.6, which corresponds to p ~ 2.2 
(Caraveo et al. 2003). Assuming that the Geminga pulsar 
moves in the interstellar medium of density n ^ 1 cm~* 
with a proper motion velocity of Vp ~120 km s^^ (Big- 
nami & Caraveo 1993), and a spin down power of Lgd — 
3.2 X 10*'' erg s~^, the termination shock radius of the pul- 
sar wind nebula can be determined as _Rs --^ 4 x lO'® cm, 
with to ~ 3 X 10® s. Taking 7^, ~ 10®, te ~ 0.5, en ~ 0.01, 
and B ~ lO"'^ G, Cheng et al. (2004b) explained the X-ray 
features of Geminga with a one-zone model. The SSC con- 
tribution to the TeV flux of Geminga is also very small. The 
total ICS luminosity of the nebula in the TeV energy range 
is estimated as ~ 3 x 10^^ erg s~^, which is lower than the 
upper limit inferred by observations. 



5 SIMULATION OF 7-RAY PULSARS IN THE 
GALAXY 

In this section we briefly describe the Monte Carlo method, 
which is used to simulate the properties of the luminosity, 
spatial evolution, and distributions of 7-ray pulsars in the 
Galaxy. The detailed Monte Carlo steps can be found in 
Cheng & Zhang 1998; Zhang et al. 2000; Fan et al. 2001). 

5.1 Monte Carlo simulation of 7-ray pulsars in 
the GalcLxy 

The basic assumptions of our Monte Caxlo simulation are 
given as follows: 

1. The birth rate of pulsars in the Galaxy is not clear, 
it is about one pulsar every 50-200 years. Here we use the 
pulsar birth rate in the Galaxy is one pulsar per century. 
The age of the Gould Belt is ~ 3 x lO'^ yrs. Again the birth 
rate of pulsars in the Gould Belt is not certain but a rate 
of 20 Myr~^ is generally used (Grenier 2000). Basically the 
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birth rates will only affect the final number of 7-ray pulsars 
but not their distributions. 

2. The Sun is inside the Gould Bolt and is located about 
200 pc towards / = 130° (Guillout et al. 1998). The Gould 
Belt has an ellipsoidal shaped ring with semi-major and mi- 
nor axes equal to 500 pc and 340 pc respectively. 

3. The initial position for each pulsar in the Galaxy is 
distributed according to the mass distribution of the Galaxy. 
We use the mass distribution suggested by Paczynski (1990) 
and Sturner & Dermer (1996). But the initial position of 
each pulsar is assumed to be born uniformly inside the Gould 
Belt. 

4. The initial magnetic field log B is assumed to be 
Gaussian, with a mean value of 12.5 and a dispersion of 0.3. 
Since the pulsars in our simulation sample are all younger 
than 10 million years, and the magnetic field docs not decay 
in 10 Myr (Bhattacharya et al. 1992), we have ignored the 
field decay for these pulsars. 

5. The initial period is chosen to be Po = 10 ms, and the 
period at time t is given by P(t) = (Po + 1.95x lO'^^S^t)^/^. 

6. The initial velocity of each pulsar is the vector 
sum of the circular rotation velocity at the birth loca- 
tion and the random velocity from the supernova explosion 
(Paczynski 1990). The circular velocity is determined by the 
Galactic gravitational potential and the Maxwellian three- 
dimensional root-mean-squarc (rrns) velocity is assumed to 
be \/3 X 100 km s~^ (Lorimer et al. 1997). Furthermore, 
the pulsar position at time t is determined following its mo- 
tion in the Galactic gravitational potential. Using the equa- 
tions given by Paczynski (1990) for the given initial veloc- 
ity, orbit integrations arc performed using the fourth-order 
Runge Kutta method with variable time steps on the vari- 
ables R,Vr,z, and Vz. The sky positions and distances to 
the simulated pulsars can then be calculated. 

7. The inclination angle a of each pulsar is chosen ran- 
domly from a uniform distribution (Biggs 1990). 

8. The 7-ray background is non-uniform over the sky, 
therefore the 7-ray threshold varies over the sky as well. In 
general the threshold will be higher in the galactic plane 
and lower at higher latitudes. It is generally accepted that 
a detectable source should have a signal at least 5a above 
the background, which is roughly equivalent to the likeli- 
hood criterion \/TS > 5(Hartman et al. 1999). Recently, 
Gonthior et al. (2002) have estimated that photon flux 
thresholds for EGRET are 1.6 x lO^^cm"^ s"^ for sources 
located at |6| < 10° and 7 x 10~*cm~^ s~^ for sources 
located at |6| > 10° respectively. In our simulation, we 
use energy flux threshold instead of photon flux threshold. 
For simplicity, we have assumed photon index ~ 2 for all 
simulated sources. We convert the photon flux thresholds 
suggested by Gonthier et al. (2002) to energy flux thresh- 
old as ~ 1.5 X 10~^°erg cm~^ s~^ for sources located at 
|6| < 10° and ~ 6.8 x 10~^^erg cm~^ s~^ for sources located 
at |6| > 10° respectively. 

5.2 Simulation results 

We carry out Monte Carlo simulation of the Galactic pulsars 
born during the past 10 Myr. We find a total of 76 7-ray ma- 
ture pulsars of ages larger than 10® yeaxs that could be de- 
tected by EGRET. Out of this simulated sample, 44 of them 
lie in the Galactic plane (|6| < 5°) and 32 lie at higher lati- 



tudes (|6| > 5°). This current result is slightly different from 
the simulations presented in Cheng et al. (2004a) because 
of the revised ago restrictions hero of taking only mature 
pulsars of > 10® years. Currently, four radio pulsars with 
age > 10® yrs are identified as 7-ray pulsars, i.e. Geminga, 
PSR B1055-52, PSR B1951-h32, PSR J0218-I-42 (Thompson 
et al. 1996; Kuiper et al. 2000). The predicted 7-ray pulsar 
numbers appear very much larger than the confirmed 7-ray 
pulsars. We should notice that first, not all 76 predicted are 
radio- loud. The radio beaming factor is roughly 0.15 (Em- 
mering and & Chevalier 1989; Biggs 1990). Taking the radio 
beaming factor into account, we predict that ~12 EGRET 
Unidentified Sources will be identified in the radio band in 
future. In fact recently there are 20 known radio pulsars lo- 
cated within the error boxes of EGRET sources and many of 
them may be identified as radio-loud 7-ray pulsars (Manch- 
ester 2004). It is very important to note that TeV 7-rays 
from pulsar wind nebulae are isotropic and independent of 
radio-loud or radio-quiet 7-ray pulsars. 

In Figure 1, we plot the distributions of the pulsar pe- 
riod (upper panel) and the magnetic field (bottom panel) of 
our simulated sample, for both sources in the the high lati- 
tudes of |6| > 5° (solid histogram), and in the low latitudes 
|6| < 5° (dashed histogram). The proper motion velocity 
(top panel) and distance (bottom panel) distributions of our 
mature pulsar sample are presented in Figure 2. The average 
velocity of our sample is found to be ~ 350 km s~^. We have 
taken this value in equation (5) to estimate the termination 
radius (cf. § 2.2). We show also in Figure 2 that the low and 
high latitude pulsar samples have different distance distri- 
butions. The average distance of the pulsars in the Galactic 
plane (|6| < 5°) is determined to be < d >~ 900 pc, while 
the average distance is < d >~ 400 pc for the high latitude 
sample (|6| > 5°). This suggests that many high latitude 
7-ray mature pulsars may lie in the Gould Belt. 

In Figure 3, we plot the 7-ray luminosity versus spin 
down power of our simulated 7-ray pulsar sample. The spin 
down power are all found to be lower than 10^® erg s~^, 
which is consistent with our expectations. Many of them in 
the Galactic plane may be the Geminga-like pulsars. If we fit 
the simulated relation of 7-ray luminosity versus spin down 
power with single power law and it is consistent with the ob- 
served correlation of the known 7-ray pulsars, L^xicV oc L^j® 
(Thompson 2001). However, from Figure 3 there probably 
exist two populations of the relations: one population fol- 
lows L-y oc Lsd; another population emerges as a branch for 
L^d > lO'^^org s"\ with oc L]^/ . Zhang et al. (2004) 
have provided a possible explanation for the nature of these 
two populations and they estimate the separation of the two 
populations should occur at 1.5 x 10^'^P^^^erg s~^. 



6 TEV SOURCE CANDIDATES OF THE 
UNIDENTIFIED EGRET SOURCES 

Assuming that mature pulsars can form faint and compact 
nebulae, and can produce TeV photons through ICS pro- 
cesses, the ICS flux from the pulsar nebulae can be calcu- 
lated by 

F.^,TeV = L^°^/i7r(fE-^, (13) 
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where L ~ Lie + issc (see § 3), and E-y is the threshold 
energy, taken to be 1 TeV. In order to calculate the TeV lu- 
minosity for each of the simulated pulsars, we have to choose 
some fixed values for the pulsar wind parameters, i.e. Se, es, 
p and fw , which are produced in non-linear processes and 
will certainly vary from pulsax to pulsax. We have assumed a 
set of pulsar wind nebula parameters for all simulated pul- 
sars, with Ce ^ 0.5 according to the energy equipartition 
assumption, es ~ 0.01, p ~ 2.2, jyj ~ 10^, and the ISM 
density n ~ 1 cm~^. The choices of these parameters are 
consistent with theoretical estimations. 

In Figure 4, we plot the GeV 7-ray photon flux from 
the pulsar outer gap versus the TeV photon flux from the 
pulsar wind nebula. There exists a correlation between the 
GeV and TeV fluxes because both quantities depend on the 
pulsar spin down power. This suggests that strong EGRET 
sources may be potential TeV sources to be detected by 
present and future TeV telescopes. 

The TeV flux distribution of the mature pulsars at high 
latitudes, \b\ > 5° (solid histogram) and on the Galactic 
plane, |6| < 5° (dashed histogram), are presented in Fig- 
ure 5. The predicted TeV fluxes of our low latitude simu- 
lated sample are all lower than 3 x 10^^^ photon cm~^ s~^, 
which is also the upper flux limit of all but two of the high 
latitude sample. The predicted TeV flux from our sample 
is lower than the previous observational constraints men- 
tioned earlier in the Introduction. However, with the rapid 
advancements of the ground-based TeV telescopes, some of 
the unidentified EGRET sources could be identified as TeV 
sources in the future. 



7 SUMMARY AND DISCUSSION 

We study the high energy radiation from mature pulsars 
with ages of ~ 10® — 10^ years in the present paper. We con- 
sider 100 MeV to GeV 7-rays generated from the magneto- 
sphere of mature 7-ray pulsars based on a new self-consistent 
outer gap model, which includes the effects of inclination an- 
gle and the average properties of the outer gap (Zhang et 
al. 2004). The relativistic wind particles from these mature 
pulsars interact with the interstellar medium, and form com- 
pact wind nebulae. The wind nebulae then produce X-ray 
emission through synchrotron radiation, and TeV photons 
from ICS by relativistic electrons on cosmic microwave back- 
ground and synchrotron seed photons. We investigate the 
nebula radiation using a one-zone model, which has been 
shown to be able to predict the high energy properties of 
known pulsar wind nebulae (Cheng et al. 2004b). We con- 
clude that in mature pulsars, the synchrotron radiation and 
the ICS processes like SSC cannot significantly contribute 
in the 100 McV to GeV energy band because of the weak 
magnetic field present {B ~ lO"'' G). However, it is possible 
that the nebulae of young pulsars in the galactic plane sur- 
rounded by supernova remnants, e.g. the Crab nebula, can 
contribute to the 100 MeV to GeV energy band through 
synchrotron and SSC processes. 

We would like to remark that pulsed TeV emission is 
expected to be produced inside the light cylinder via inverse 
Compton scattering between the soft photons and the rel- 
ativistic charged particles inside the outer gap. In the new 
outer gap models (Zhang et al. 2004), the soft photon density 



is just enough to convert one out of 10® gamma-ray photons 
into pairs. The inverse Compton scattering efficiency is ex- 
pected to be lower than 10~® of the outer gap power. There- 
fore, it is not surprised that none of known EGRET pulsars 
has been detected with pulsed TeV photons. On the other 
hand, TeV photons associated pulsaxs should be produced 
by PWN via the inverse Compton scattering and therefore 
should be unpulsed. The Crab nebula has been confirmed 
with TeV detection. PSR1706-44 is considered as another 
confirmed TeV source (Weekes 2004) but the recent HESS 
results challenge this claim. Geminga is a very nearby known 
EGRET pulsar but yet no positive detection on TeV. For the 
former pulsar, we have speculated that the emission from 
PWN may be variable with a time scale of years and the lat- 
ter may attribute to the fact that the magnetic field strength 
in PWN is so weak that the TeV fiux produced by SSC is 
also very weak. In short, there are some evidences that TeV 
photons can be emitted from some known EGRET pulsar 
PWNe. Therefore we beheve that if the EGRET unidenti- 
fied sources are indeed pulsars, some of their PWNe should 
emit detectable TeV photon fluxes. 

Finally, we calculate the TeV fluxes from wind nebulae 
of these mature pulsars through inverse Compton processes. 
The predicted TeV fluxes axe consistent with the recent ob- 
servational constraints obtained at the Whipple Observa- 
tory. In addition, our results predict that the unidentified 
EGRET sources, especially the strong sources, can be po- 
tential TeV sources to be be detected by future ground-based 
TeV telescopes. Several third generation ground-based TeV 
telescopes are under construction or recently finished, e.g., 
MAGIC (Lorenz 1999; 2004), HESS (Hofmanu 1999; Hin- 
ton 2004), VERITAS (Weekes et al. 2002; Krennrich et al. 
2004), and CANGAROO-III (Matsubara 1997; Kubo et al. 
2004). The high sensitivity of these new TeV telescopes of- 
fer the possibility of investigating the unidentified EGRET 
sources. For example, HESS is an imaging telescope array 
system with an array of four imaging atmospheric Cherenkov 
telescopes. Now it is fully operational, and its expected sensi- 
tivity at 1 TeV is about 2 x 10~^^ photon cm~^ s~^, then we 
also expect HESS could possibly detect TeV photons from 
about 15 unidentified EGRET sources if they are mature 
pulsars. 

We thank the anonymous referee for his very useful 
comments and Prof. P.K. MacKeown for his critical read- 
ing. This work is supported by a RGC grant of the Hong 
Kong Government. 
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Figure 1. The distribution of the period (upper panel) and sur- 
face magnetic field (bottom panel) of the simulated 7-ray pulsars 
which could be detected by EGRET. The distributions of the pul- 
sars in the high latitude (|6| > 5°) (solid), and in the Galactic 
disk (|6| < 5°) (dashed) are shown. 



Figure 2. The distribution of the proper motion velocity (top) 
and the distance (bottom) of the simulated 7-ray pulsars. The 
distributions of the pulsars in the high latitude {\b\ > 5°) (solid), 
and in the Galactic disk (|6| < 5°) (dashed) are shown. 



© 0000 RAS, MNRAS 000, 000-000 



Possible TeV Source Candidates Among The Unidentified EGRET Sources 11 



"> 33.5 



<? 33.0 



32.5 




32.0 32.5 33.0 33.5 34.0 34.5 35.0 35.5 36.0 

Log (erg s ') 



E 6 



"T 1 1 ' 1 1 1 ' 1 1 r 



-14.5 -14.0 -13.5 -13.0 -12.5 -12.0 -11.5 -11.0 -10.5 

Log L (photon cm"^ s"') 



Figure 3. The GeV gamma-ray luminosity versus the spin down 
power for the simulated 7-ray pulsars for the high latitude |6| > 5° 
(solid) and the Galactic latitude |6| < 5° (circle) sample. 



Figure 5. The distribution of the TcV flux from the wind nebu- 
lae of the simulated 7-ray pulsars which could be the unidentified 
EGRET sources. The distributions of the pulsars in the high lati- 
tude (|6| > 5°) (sohd), and in the Galactic disk (|6| < 5°) (dashed) 
are shown. 
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Figure 4. The GeV gamma-ray flux versus the TeV flux for the 
simulated 7-ray pulsars for the high latitude |6| > 5° (solid) and 
the Galactic latitude |6| < 5° (circle) sample. 
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